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ABSTRACT 

An  analysis  is  presented  which  predicts  the  properties  of  an  arbitrarily  thick  turbulent  boundary 
layer  in  axial  fiow  past  a  long  cylinder.  The  study  makes  use  of  a  modified  form  of  the  turbulent 
law-of-the-wall,  which  properly  accounts  foi  transverse  curvature  effects.  Using  this  law,  the 
theory  which  follows  is  then  an  exact  solution  to  the  axisymmetric  equations  of  continuity  and 
momentum  in  incompressible  flow. 

Numerical  results  are  given  to  show  the  effect  of  curvature  on  the  various  boundary  layer  char¬ 
acteristics.  Skin  friction  and  drag  coefficients  can  be  increased  greatly  with  increasing  curvature 
while  boundary  layer  thickness  isdecreased.  When  defined  in  their  axisymmetric  form,  the  displace¬ 
ment  and  momentum  thickness  are  both  decreased  b^  curvature.  The  velocity  profile  is  flattened 
greatly  and  the  shape  factor  H  -i*/9  approaches  unity  at  large  curvature.  The  failure  oi  earlier 
power-law  theories  to  make  accurate  predictions  is  shown  to  be  due  to  their  inadequate  handling 
of  the  strong  profile  shape  changes.  Finally,  the  earlier  concept  that  the  curvature  effect  on  skin 
friction  could  be  correlated  by  the  ratio  (d/d)  is  shown  to  be  invalid. 
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NOMENCLATURE 


a  - 
B,  k 

Cf 

CD 

d 

f 

H 

L 

r,  x 

Rx »  Rl  = 

u,  V 

U 

v*  = 

y 

Y 
i 

i*,  9  = 

P 

M  = 

V  = 

r  = 

X 


cylinder  radius 

law-of-the-wall  constants,  Eq.  (2) 

local  friction  coefficient  =  2rw/pU2 

overall  drag  coefficient  =  2{Dx*ag)/pU22raL 

cylinder  diameter  =  2,i 

iaw-of-the-wall  function,  Eq.  (12) 

shape  factor  =  i*/0 

cylinder  length 

radial  and  axial  coordinates 

local  and  overall  Reynolds  numbers 

axial  and  radial  velocity  components 

free  si  ream  velocity 

wall  friction  velocity  =  ^rw/p 

wall  coordinate  =  [r-a] 

modified  law-of-the-wall  variable,  Eq.  (10) 

boundary  layer  thickness,  Fig.  1 

displacement  and  momentum  thickness,  Eqs.  (24),  (25) 

fluid  density 

fluid  absolute  viscosity 

fluid  kinematic  viscosity  =  m/p 

boundary  layer  shear  stress 

dimensionless  friction  variable  =  U/v*  =  ^2/Cf 


Subscript: 

w  =  at  the  wall 
Superscript: 

+  =  iaw-of-the-wall  variable,  Eqs.  (1),  (7) 


THE  AXEYMMETRJC  TURBULENT  BOUNDARY  LAYER 
ON  AN  EXTREMELY  LONG  CYLINDER 


INTRODUCTION 

The  prediction  of  boundary  layer  properties  in  axial  How  past  long  bodies 
of  revolution  is  an  important  practical  problem.  Such  bodies  are  common  in 
both  aerodynamics  and  hydrodynamics,  and  one  wishes  to  know  whether  the 
transverse  curvature  has  any  important  effect  in  estimating,  say,  skin  friction 
or  displacement  thickness.  If  the  flow  is  turbulent,  the  prediction  of  wall  pres¬ 
sure  fluctuations,  along  with  the  coupled  dynamic  response  that  such  flow  noise 
might  cause  to  the  body,  is  important.  As  a  first  step  in  such  problems,  it  is 
the  purpose  of  this  report  to  analyze  the  axial  flow  turbulent  boundary  layer 
past  an  extremely  long  circular  cylinder.  The  analysis  is  based  upon  a  cylin¬ 
drical  version  of  the  turbulent  law -of-the -wall ,  first  deduced  by  G.  N.  V.  Rao.13 
The  present  study  is  limited  to  incompressible  flow  with  constant  properties. 


LITERATURE  REVIEW 

First  let  us  remark  that  the  case  of  laminar  flow  past  a  long  cylinder  has 
been  dealt  with  very  successfully.  The  laminar  boundary  layer  equations 
achieve  closure  without  any  extra  empiricism  needed,  and  accurate  analytical 
solutions  have  been  given  by  Probstein  and  Elliott1  and  byGlauert  and  Lighthill.2 
Other  laminar  flow  theories,  mostly  introductory  in  nature,  are  given  by  Seban 
and  Bond,3  Kelly,4  and  Cooper  and  Tulin.5  Very  strong  curvature  effects 
occur  in  laminar  flow.  For  example,  at  x/d  =  100,  the  local  skin  friction  is 
increased  nearly  100  percent  over  a  flat  plate  with  the  equivalent  Reynolds 
number.  A  theory  for  supersonic  laminar  cylinder  flow  was  given  by  Mark.® 


As  might  be  expected,  the  bulk  of  experimental  data  occur  under  turbulent 
flow  conditions,  which  is  usually  the  case  also  in  practical  applications.  Prob¬ 
ably  the  earliest  data  on  axial  cylinder  flows  are  due  to  Kempf ,  ?  and  subsequent 
experiments  have  been  performed  by  Telfer,8  Hughes,9  Richmond,10  Yu, 11 
Yasuhara, 12  and  Rao. 13  The  data  of  Richmond10  are  the  most  extensive  and 
encompass  both  incompressible  and  hypersonic  turbulent  flows.  The  experi¬ 
ment  of  Yasuhara12  included  laminar  flow  conditions  and  showed  transition  to 
turbulence  at  a  Reynolds  number  of  R*  =  1.5  x  10®,  not  much  different  from 
typical  smooth  flat  plate  transition  points. 


The  earliest  theory  for  turbulent  cylindrical  flow  is  due  to  Millikan, 14  who 
established  the  proper  form  of  the  Karman  integral  relation  for  axisymmetric 
flow.  Millikan's  calculations,  which  use  the  familiar  1/7  power  ’aw  for  the 
velocity  profile,  do  not  agree  with  experiment,  because  the  1/7  power  is  not  at 
all  accurate  when  curvature  effects  exist.  Further  theories  using  power-law 
profiles  were  given  by  Landweber,  15  Eckert,  16  Karhan,  17  and  Sakiadis.  16 
These  results  do  not  agree  with  experiment  either  and  generally  predict  curva¬ 
ture  effects  that  are  too  small  by  an  order  of  magnitude. 

A  rather  different  approach  was  taken  by  Ginevskii  and  Solodkin,  1°  who 
adapted  Prandtl's  mixing  length  theory  to  the  thick  cylindrical  boundary  layer. 
Their  theory,  which  is  extremely  complex  algebraically,  includes  both  concave 
and  convex  surfaces  and  pressure  gradient  and  flow  separation  effects.  Unfor¬ 
tunately,  their  mixing  length  parameters  were  taken  from  flat  plate  data  and 
thus  are  not  accurate  under  strong  curvature  conditions. 


THE  LAW-OF-THE-WALL  FOR  A  CYLINDER 

As  is  well  known,  the  equivalent  shear  stress  in  turbulent  flow  is  not  simply 
related  to  viscosity  and  velocity  gradient  but  includes  a  large  fluctuating  inertia 
term,  the  physics  of  which  is  not  well  defined  even  to  the  present  day.  How¬ 
ever,  in  two-dimensional  flow,  under  flat  plate  (zero  pressure  gradient)  condi¬ 
tions,  it  is  a  fortunate  circumstance  that  the  local  velocity  u(x,y)  may  be 
related  with  great  accuracy  to  local  wall  variables:  u  --  fcn(y,  p ,  p ,  rw) .  In 
dimensionless  form  this  becomes  the  celebrated  law -of-the -wall: 

u+  =  u/v*  =  fcn(y+) ,  (1) 

where  y+=  Pv* y/p  and  v*  =  yJT^/p  •  Very  close  to  the  wall,  viscous  shear  is 
dominant  and  a  sublayer  occurs,  where  u+  »  y+  .  Away  from  the  wall,  turbu¬ 
lent  shear  predominates,  resulting  in  a  logarithmic  region: 

u+  =  £  I"  (y+)  +  B  ,  (2) 


where  k  and  B  are  empirical  constants  that  hold  accurately  for  all  flat  plate 
or  modest  pressure  gradient  flows.  Coles120  suggests  (k  =  0.41,  B  =  4.9), 
while  Spalding21  prefers  the  nearly  equivalent  values  (k  =  0.40,  B  =  5.5).  The 
difference  is  not  important. 


Finally s  very  far  from  the  wall,  a  wakelike  flow'  (Coles30)  occurs  which  is 
independent  of  viscosity.  We  shall  not  consider  this  wake  in  the  present  report, 
because  numerous  theories  have  shown  that  the  outer  wake ,  while  very  inter¬ 
esting  physically,  has  a  negligible  numerical  effect  on  boundary  layer  computa¬ 
tions  such  as  skin  friction  arid  displacement  thickness.  / 

The  law-of-the-wall  provides  the  necessary  closure  relation  for  the  turbu¬ 
lent  boundary  layer,  relating  wall  shear  to  u  and  y.  No  further  assumptions 
are  needed  for  incompressible  flow,  and  Eq.  (2)  may  be  combined  with  the  con¬ 
tinuity  and  momentum  relations  to  yield  very  accurate  and  simple  flat  plate 
turbulent  calculations,  as  in  the  theory  of  Brand  and  Per  sen.  22 

For  cylinder  flow,  the  entire  concept  must  be  reevaluated.  Here  the  cyl¬ 
inder  radius  becomes  an  additional  parameter,  and  we  may  postulate  that 
u  =  fcn(r , a,  p ,  m  ,  rw),  leading  to  any  of  several  possible  dimensionless  forms; 
for  example, 

*4* 

u  -  fen  (a  ,  r/a)  (Rao13) 

or  +  +  .  (3) 

u  =  fcn(y  ,  y/a)  (Richmond10), 

where  y  is  the  wall  coordinate  such  that  r  =  y  +  a .  See  Fig.  1  for  a  definition 
sketch  of  terms.  The  formulation  of  Richmond  above,  the  earlier  of  the  two  by 
ten  years,  was  based  upon  the  so-called  "streamline  hypothesis"  of  Coles. 20 
Coles  was  the  first  to  show  that,  in  two-dimensional  flow,  lines  of  constant  u+ 


Fig.  i.  Definition  Sketch  for  Axial  Cylinder  Flow 
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are  streamlines  of  the  flow.  Richmond*0  hypothesized  that  the  same  would 
be  true  in  cylindrical  flow,  where  the  stream  function  varies  as  (r2  -  a2). 

If  this  were  true,  uT  would  then  be  a  function  of  only  a  single  parameter 
[(a  +  y)2  -  a2] ,  instead  of  a  and  y  separately.  Thus,  Richmond  proposed 

that 

u+ =  fcnjv*(r2  -  a2)/2avj  =  fcn|y+(l  +  y/2a)J  .  (4) 

This  relation  is  also  used  to  estimate  v*  from  the  data.  Upon  plotting  the 
results  from  various  experiments  in  th*«  manner  —  see  Fig.  7  of  Yasuhara*2 — 
the  curves  do  indeed  correlate  very  near  the  wall.  Away  from  the  wall, 
thougu,  the  curves  begin  to  drop  off  significantly  from  the  log  )aw,  Eq.  (2). 
Furthermore,  the  deviations  are  a  rather  erratic  function  of  the  relative  bound¬ 
ary  layer  thickness  (5/a).  These  facte  tended  to  discourage  further  analytical 
work  along  lrrw  -of -the -wall  lines,  until,  in  1967,  Rao*3  pointed  out  the  flaw  in 
Richmond's  approach. 

The  critical  point  made  by  Rao  is  that  Richmond's  proposal,  Eq.  (4), 
fa'ls  to  describe  the  sublayer  properly.  Equation  (4)  leads  to  the  sublayer  pro¬ 
portionality 

u+  =  y+  (1  +  y/2a)  (Richmond*0) .  (5) 

However,  in  the  sublayer-,  it  was  shown  by  Rao  that  the  true  relation  is 
quite  different.  Inertia  is  negligible  in  the  sublayer,  and  hence  the  boundary- 
layer  momentum  equation  for  zero  pressure  gradient  reduces  to 


Integrating  this,  we  obtain  the  proper  sublayer  relation 

u+  =  a+  In  (r/a) ,  (7) 

where  a+  =  v*aA  .  ^he  power  series  approximation  of  the  logarithm  for  small 
y  leads  to  the  relation 

u+  =  y+  (1  -  y/2a)  (Rao*3) .  (8) 
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The  difference  in  sign  of  the  parenthesis  terms  in  Eqs.  (5)  and  (8)  is  profound 
when  the  sublayer  is  thick,  which  of  course  it  is  whenever  curvature  is  impor¬ 
tant.  T!  is,  Richmond's  expression  fails  just  when  we  need  it  most:  in  esti¬ 
mating  the  wall  shear.  This  is  shown  by  solving  Eqs.  (5)  and  (7)  for  the  wall 
shear: 


Richmond:  r  =  Mu/y  {1  +  y/2a) 

W  (9) 

Rao  (exact):  r  =  fi u/a  in  (1  +  y/a) . 

w 

Thus,  Richmond's  expression  underestimates  the  wall  shear  by  the  factor 
[(y/a)  (1  +y/2a)/  In  (1  +  y/a)]  .  The  difference  is  negligible  if  (y/a)  is  very 
small  (thin  boundary  layer),  but  in  Richmond's  thick-layer  experiments,  the 
measurement  closest  to  the  wall  was  at  about  y/a  =  0. c  The  error  involved 
is  the  factor  [0.5  (1.25)/in  (1.5)]  =  1.54;  the  true  shear  was  actually  54  per¬ 
cent  larger  than  Richmond's  tabulated  values. 


From  his  sublayer  estimate,  Eq.  (7),  Rao13  then  hypothesized  that  the 
cylindrical  law -of -the -wall  might  actually  satisfy  the  usual  two-dimensional 
formulas,  if  y+  is  replaced  by  the  new  sublayer  variable  Y  =  a+  in  (r/a). 
For  example,  the  logarithmic  layer  on  a  cylinder,  using  Spalding's  constants, 
would  be  given  by 


u+  =  2.5  In  ( Y)  +  5.5,  Y  =  a+  Jn  (r/a) . 


(10) 


Figure  2  shows  all  available  data  on  thick  cylindrical  turbulent  boundary 
layers,  plotted  in  the  suggested  form  u+  versus  Y.  It  is  seen  that  Rao's  hy¬ 
pothesis  is  indeed  the  correct  one.  The  agreement  is  good  enough,  considering 
the  wide  variety  of  relative  curvatures,  so  that  we  may  regard  Fig.  2  as 
the  correct  closure  relation  needed  to  complete  the  analysis  of  cylindrical 
flows. 
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THE  THICK  TURBULENT  CYLINDRICAL  BOUNDARY  LAYER 
SKIN  FRICTION  AND  DRAG 

The  modified  law-of-the-wall  from  Fig,  2  is  the  missing  link  needed  to 
complete  the  theory  of  turbulent  flow  past  an  extremely  long  cylinder.  For  in¬ 
compressible  axisymmetric  flow  with  zero  pressure  gradient,  the  two  basic 
relations  are  the  equations  of  continuity  and  momentum: 

Continuity:  &&  +  =  0 

a*  ^ 

(id 

Momentum:  ru^  +  rv—  =  —  —  (rr) . 

ax  ay  p  ay 

If,  now,  the  law-of-the-wall  is  given  by  the  functional  relation 

u(x,y)  =  v*(x)  f(Y) ,  (13) 

then  Eqs.  (11)  and  (12)  may  be  solved  directly  for  the  wall  shear.  The  velocity 
v  may  be  eliminated  through  continuity: 

's  y  „  o  /W 

rv^-^rrudy^-rf— .  (13) 


Substitution  of  u  and  v  from  Eqs.  (12)  and  (13)  gives  the  following  differential 
equation  in  the  wall  shear  stress  alone: 


which  we  may  integrate  from  the  wall  to  the  outer  edge  of  the  boundary  layer: 


r  =  -  f  f2  (r/a)2  dY . 


av*  c 
'dx  j 


This  relation  is  essentially  identical  to  the  flat  plate  skin  friction  expression 
derived  by  Brand  and  Persen,^2  except  for  the  additional  "curvature"  term 
(r/a)2  =  exp  [2Y/a+] .  Clearly  the  wall  shear  is  increased  by  curvature,  the 
effect  being  negligible  if  r  =  a  (thin  boundary  layer). 
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Let  us  now  define  a  dimensionless  skin  friction  variable  related  to  the 
usual  coefficient  Cf .  Let 


X  =  U/v*  =  j27Cf, 


(16) 


where  Cj  =  2rw/pU2  .  Substitution  into  Eq.  (15)  gives  a  relation  between  the 
local  Reynolds  number  and  the  skin  friction  coefficient,  which  is  the  central 
result  of  this  report: 


Ux/*>  =  R  = 


A 

{ 


Y(X)  2 

G{X)  dX  ,  G  =  f  r  exp  [2Y/a  ]  dY .  (17) 


J 


It  merely  remains  to  carry  cut  the  double  integration  indicated  by  Eq.  (17),  by 
using  some  suitably  accurate  representation  of  the  modified  law-of-the-wali, 
u+=  f(Y) .  The  logarithmic  law,  Eq.  (10),  is  reasonably  accurate  at  large  Rx , 
provided  care  is  taken  to  eliminate  the  unwanted  effects  of  the  singularity  at 
the  wall.  However,  the  sublayer  becomes  increasingly  important  with  increased 
curvature.  Hence,  for  maximum  accuracy,  we  choose  to  work  with  the  more 
cumbersome  relation  of  Spalding,2!  which  fits  all  three  wall  regions  —  sublayer, 
buffer  layer,  and  log  layer— very  accurately: 

Y  =  u+  +  0.1108  [eZ  -  1  -  z  -  z2/2  -  z3/6  -  z4/24]  ,  (18) 

where  z  =  ku+  =  0.4  u+. 

The  constant  0.1108  corresponds  to  the  choice  B  =  5.5  in  Eq.  (2).  Equa¬ 
tion  (18)  is  plotted  for  comparison  in  Fig.  2,  and  the  agreement  in  all  regions  is 
seen  to  be  good.  Note  that  we  require  the  modified  variable  Y  =  a+  in  (1  +  y/a), 
not  y+.  Since  Eq.  (18)  is  a  bulky  and  implicit  form  of  f(Y),  its  use  in  Eq.  (17) 
makes  numerical  integration  a  necessity.  The  author's  numerical  results, 
which  were  obtained  by  use  of  a  fourth-order  Runge-Kutta  procedure,  are 
shown  in  Fig.  3  as  a  plot  of  Of  =  2/X2  versus  Rx ,  with  the  position  (x/d)  on 
the  cylinder  as  a  parameter.  The  original  parameter,  a+ ,  is  eliminated 
through  the  relation  x/d  =  Rx/2a+X .  Note  that  realistic  values  of  curvature 
(x/d  =  105)  may  increase  the  local  skin  friction  by  a  factor  of  ten  or  more; 
this  result  is  made  possible  by  the  profound  changes  in  the  velocity  profile 
implied  by  the  modified  law -of-the -wall,  Eq.  (12). 
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Fig.  3.  Turbulent  Skin  Friction  on  ■  Long  Cylinder, 
fro«  Eq».  (17)  end  (18) 

The  total  drag  on  the  cylinder  is  evaluated  by  integrating  the  shear  stress, 
and  the  drag  coefficient  Cq  is  defined  with  respect  to  the  wetted  area: 


Drag  = 


L, 

S 


r  2ira  dx  , 
w 


and 


(19) 


CD  "  ..,2 


2  Drag 


PU  2raL 


Drag  calculations  were  performed  in  two  ways  and  were  checked  against  each 
other  to  four  significant  figures.  First,  from  the  definition  of  Cf,  Eq.  (19) 
may  be  rewritten  as 

1 

CD  =  |  Cf  d(x/L)  .  (20) 

0 


9 


Second,  if  the  leading  edge  is  a  uniform  flow  u  =  U  ,  the  drag  may  be  calcu¬ 
lated  by  a  momentum  integral  across  the  trailing  edge: 


CD  »  (2/L) 


[u/U{l  -  u/U)  (l  +  y/a>]x  =  . 


where  the  velocity  profile  is  found  from  the  modified  law-of-the-wall,  Eq.  (18), 
and  the  known  value  of  Cf  at  the  trailing  edge.  The  two  approaches  should 
give  identical  results.  The  numerical  drag  results  are  plotted  in  Fig.  4  as  a 


\jo4 
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Drag  Coefficient  of  a  Long  Cylinder  Compared  with  That  of 
a  Flat  Plate  at  the  Same  Reynolds  Number 


ratio  of  Cq  on  the  cylinder  to  the  equivalent,  value  for  a  flat  plate,  which  is 
correlated  accurately  by  the  well  known  formula  (see,  for  example, 
Schlichting,  23  £q.  (21.16)). 

CD(flat  plate)  =  0 . 455/tlog^Rj2  * 58 .  (22) 

Figure  4  shows  that  the  drag  of  a  long  cylinder  may  be  up  to  five  times  the 
equivalent  plate  drag.  Drag  ratios  higher  than  five  are  probably  not  realistic, 
because  the  viscous  sublayer  then  extends  almost  entirely  across  the  boundary 
layer  and  raises  stability  questions  about  the  assumption  of  turbulence. 


VELOCITY  PROFILES  AND  INTEGRAL 
THICKNESSES 

Although  the  drag  and  skin  friction  of  the  cylinder  are  the  primary  results 
of  this  report,  certain  applications  —  e.g. ,  wall  pressure  fluctuations  —  are 
related  to  the  shape  and  thickness  of  the  velocity  profiles. 


Once  Cf  =  2/A2  is  known,  the  boundary  layer  thickness  is  specified  by  the 
modified  law -of-the -wall:  A  =  f[Y(5)].  That  is,  X  determined  Y(i)  and  then 
6  is  calculated  from 


5/a  =  exp  [Y(5)/a+]  -  1  . 


(23) 


Similarly,  the  law-of-the-wall  determines  r  =  u+A  »  so  that  the  profile 
shape  is  a  function  of  two  parameters:  (Cf,u)  or,  better,  (Rx ,  x/d) .  Cur¬ 
vature  causes  the  profile  to  be  much  flatter,  and  typical  shapes  (for  Rx  =  10  ^) 
are  shown  in  Fig.  5.  Also  shown  are  values  of  the  exponent  n  that  make  the 
curves  fit  approximately  to  the  familiar  power  law  (u/U)  =  (y/5)n.  We  see 
that  the  curvature  changes  the  power  in  this  case  from  about  l/llth  to  about 
l/24th.  This  is  one  reason  for  the  failure  of  the  earlier  power-law  theories^-18 
to  make  accurate  predictions;  another  reason  is  simply  that  the  power  law  is 
not  really  a  good  fit  anyway,  especially  in  the  vital  sublayer  region  which  cur¬ 
vature  emphasize?. 


The  momentum  and  displacement  thicknesses  can  also  be  calculated  from 
the  law-of-the-wall.  Kelly^  was  the  first  to  point  out  that  the  definition  of 
these  thicknesses  should  change  to  account  for  the  cylindrical  geometry.  Thus, 
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the  displacement  thickness  on  a  cylinder  is  given  by  the  formula 

a+3 


1 


(6*  +  a)2  -  a2  =  |  [l  -  u/u]  2r  dr  , 


(24) 


and  the  momentum  thickness  is  given  by 

a+i 


1 


(0  +  a)2  -  a2  =  j  (u/U)  [l  -  u/u]  2r  dr  . 


(25) 


No  other  definition  is  proper  because  it  will  not  properly  model  the  cylindrical 
mass  flow  relations  along  the  cylinder.  Values  of  6  and  5*  calculated  from 
these  relations  are  shown  in  Fig.  6  for  various  Reynolds  numbers  and  curvatures. 


Fig.  6.  Boundary  Layer  and  Displacement  Thicknesses, 
from  Eqs .  (23)  and  (24) 


We  see  that  curvature  reduces  both  6  and  6*  by  approximately  the  same 
amount,  down  to  no  more  than  a  factor  of  about  ten. 

The  momentum  thickness  calculated  from  Eq.  (25)  is  only  slightly  smaller 
than  S* ,  as  would  be  expected  from  the  profile  flattening  effect  shown  in 
Fig.  5,  Figure  7  shows  values  of  the  so-called  "shape  factor"  H  =  5 */8  for 
various  Reynolds  numbers  and  curvatures.  The  shape  factor  drops  to  nearly 
unity  as  (x/d)  increases,  indicating  a  nearly  flat  (slug  flow)  profile.  There 
is  substantial  effect  even  at  the  very  small  curvature  of  x/d  =  iO  .  The  some¬ 
what  erratic  curve  for  x/d  =  10, 000  reflects  that  in  this  case  the  viscous  sub¬ 
layer  extends  nearly  across  the  entire  boundary  layer.  Except  for  a  few  cases 
near  (x/d  =  0) ,  we  may  state  that  curvature  decreases  6  just  as  it  does  6*  . 
This  seems  to  be  in  conflict  with  the  accompanying  drag  increase  until  we  re¬ 
member  that  6  ,  as  defined  by  Eq.  (25),  is  not  a  measure  of  drag,  whereas  in 
two-dimensional  flow  the  flat,  plate  drag  coefficient  Cd  =  20/L  . 


»« 


Fig.  7.  Theoretical  Valuea  of  the  Turbulent  Flow  Shape  Factor  on 
Extremely  Long  Cylinder*,  from  Eqa.  (24)  and  (25) 
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for  comparison  and  are  seen  to  be  in  excellent  agreement  with  the  present 
theory.  The  dotted  curve  shows  Richmond's  original  suggested  correlation, 
which  is  now  cleaxly  inadequate  in  view  of  the  data  reevaluation  seen  in  the  light 
of  the  modified  law-of-the-wall.  Since  Fig.  8  requires  two  parameters  to  spec¬ 
ify  a  point,  it  is  now  no  better  in  principle  than  the  presently  suggested  Fig.  3. 
In  practice,  Fig.  8  is  much  less  handy  than  Fig.  3,  because  the  calculation  of 
$  from  Eq.  (25)  is  a  cumbersome  matter,  whereas  the  estimation  of  Rx  and 
x/d  is  relatively  simple. 


CONCLUSIONS 

It  has  been  shown  that  the  modified  cylindrical  law-of-the-wall,  Fig.  2, 
leads  to  a  complete  and  rational  analysis  of  the  turbulent  boundary  layer  in 
axial  flow  past  long  cylinders.  The  present  theory  is  straightforward,  leading 
to  Eq.  (17),  and  the  numerical  results  are  in  excellent  agreement  with  available 
experimental  data.  The  theory  seems  to  have  no  apparent  limitations  on  cyl¬ 
inder  size  or  Reynolds  number,  as  long  as  the  flow  is  turbulent.  At  low  Rey¬ 
nolds  numbers  and  large  curvatures,  the  velocity  profiles  are  mostly  viscous 
sublayer,  suggesting  that  in  such  cases  the  flow'  might  actually  be  laminar. 
Tables  and  charts  of  important  boundary  layer  parameters  are  easily  con¬ 
structed  from  the  theory,  and  several  charts  are  included  here  for  the  reader '3 
interest.  Also,  the  approach  could  easily  be  extended  to  pressure  gradient 
flows  if  warranted . 
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